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ABSTRACT 

We study 87 Extremely Red Objects (EROs), selected both to have color redder than R— [3.6] = 4.0, 
and to have confirmed spectroscopic redshifts. Together, these two constraints result in this sample 
populating a fairly narrow redshift range at 0.76 < z < 1.42. The key new ingredient included here is 
deep Spitzer Space Telescope InfraRed Array Camera (IRAC) data. Based on [3.6] — [8.0] color, we 
demonstrate that it is possible to classify EROs into early-type, dusty starburst, or power-law (AGN) 
types. We present ultraviolet to mid-infrared spectral energy distributions (SEDs) and Advanced 
Camera for Surveys (ACS) images, both of which support our simple IRAC color classification. 
Subject headings: galaxies: evolution — galaxies: high-redshift — infrared: galaxies — galaxies: ellip- 
tical and lenticular, cD — galaxies: starburst 



1. INTRODUCTION 

First discovered in the late 1980s (Elston et al."1988^, 
extremely red objects (EROs) are defined by their very 
red optical/near-infrared colors. It has been known for 
some time that the redness of their color constrains these 
galaxies to be either early-type galaxies, starburst galax- 
ies reddened by dust, or AGN (or a combination of these 
three classes). However, until recently, with observations 
limited to i^T-band (or shorter wavelengths), it has proven 
extre mely challen ging to accurately classify the EROs by 
type (|Mannucci e t al. 200^), even in combina tion with 
I high-resolution HST imaging (Moustakas et al. 2004). 

Although EROs appear to consist of a heterogeneous 
mix of galaxy classes, the emerging paradigm is that 
they may well be the high redshift counterparts and pro- 
genitors of local massive E and SO galaxies. Their re- 
liable classification and study, especially at intermedi- 
ate redshift {z ^ 1) can provide crucial constraints on 
the evolution of massive, starburst, dusty and/or ultra- 
luminous infrared galaxies (ULIRGs) known to exist at 
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higher redshift e.g BzKs fPaddi et ai;'2QQ4 !), BX/BMs 
(iReddv et al.M 2005D. Distant Red Galaxies (jFranx et all 
2003; Papovich^TaUlinni Eonselice 2003, DRGs), sub- 
millimeter and IR-Lu minous Lyman Break Galaxies 
(Huang et al. 2005; Rigopoulou' 2006, ILLBGs). 

In the restframe near-IR, old stellar populations show 
a turndown at wavelengths longer than the restframe 
1.6/im 'bump', while dusty starburst populations show 
emission from small hot dust grains. AGN-dominated 
sources display a_ powe r-law spectral energy distribution. 
In Wi lson et al.l (|2004l) we showed how data from Spitzer 
could begin to help distinguish among different ERG 
populations. In this paper, we extend our ERG study to 
take advantage of the rich panchromatic dataset avail- 
able from the All-wavelength Extended Groth strip In- 
ternational Survey (AEGIS). All magnitudes used in this 
letter are AB, unless otherwise specified. 

2. THE AEGIS DATASET 

Th e Spitzer IRAG (3.6, 4.5, 5.7, 8.0/im, iFazio et al.1 
I2004D component of the Extended Groth Strip (EGS) 
survey spans an area of 120 x 10 arcmin (Huang et al. 
2006, in prep: see also Huang et al. 2004; Barmbv et al.l 
120061 IHuan g 2006). In conducting this ERG study, we 
also utilized DEEP2 spectroscopy (Davis et al. 2006, in 
prep), dnid u'q' (Ashby et al. 2006, in prep), GFHT BRI 
(ICoil et a1.! ^4). AGS V(F 606W) and I{F8UW), deep 
Subaru R (MivazakTeEaDr2002, 27.0 AB, 5cr), K (Gon- 
selice et al. 2006, in prep), and Spitzer Multiband Imag- 
ing Photometer (MIPS, Rieke et al.^2004) 24/ i m ima ging. 
Further details may be found in lDavis et"al] (|2006D . 

3. IRAC-SELECTED GALAXIES AND THEIR 
COLOR-COLOR DISTRIBUTION 

The EGS field contains ~ 45000 galaxies detected at 
3.6/im (23.9 AB, 5a). Figure d shows an — [3.6] ver- 
sus [3.6] — [8.0] color-color diagram. The black points 
show the ~ 13000 3.6/im-selected galaxies with both 
good quality Subaru R and 8.0/im photometry. 

The colored tracks on Figure Q show the location in 
color-color space as a function of redshift (0 < z < 4) 
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for ei^ht non-evolving empirical tem plates (four com- 
mon ESimanI^I^]^edma3ill^ templates 
[E, Sbc, Scd, Im] empirically extended to 10/im using 
ISO data (Huang et al., 2006, in prep), a dusty star- 
burst template [M82], a dusty starburst/ULIRG tem- 
plate [Arp220], and two AGN templates [NGC 1068 and 
NGC 5506]). While we fully expect these templates to 
become increasingly inaccurate at high redshift, they do 
serve to provide simple insight onto the likely nature and 
redshift distribution of galaxies within this color-color di- 
agram. 

A distinctive swath of galaxies is clearly apparent curv- 
ing from the left to the top of Figure (the plume in 
the lower left corner is caused by stellar contamination). 
Notice, especially, the excellent agreement between this 
swath of galaxies observed in color-color space and the 
CWW E (blue) track to z ^ 0.7, when the template 
noticeably begins to diverge from the "bluer" data. We 
note that the predictability of the [R] — [3.6] color-redshift 
relation for early-type galaxies can be utilized as an ef- 
fective technique for detecting high redshift clusters of 
galaxies e.g., at 2; < 1.4 in the 4d eg:^ Spitzer First Look 
Survey Field (|Wilson et al.ll2005l: Muzzin et al, 2006, in 
prep). The SpARCS collaboration-^^ is currently utiliz- 
ing an even redder [z'] — [3.6] color to detect and study 
cluster s to z = 2 in the 50deg^ Spitzer SWIRE Legacy 
Fields (|Wilson et al.i r2006). 

We define an ERO to be a galaxy redder in color than 
R — [3 .6] = 4.0. This is the same criterion used in 
iWilson et all (2004), and is a very similar selection cri- 
teri on to the traditional Vega R — K > 5.0 requirement 
f see IWilson et al. 2004 for a discussion). 

From Figure^ we might expect that only high redshift 
{z > 0.8) early- type, {z > 1.0) dusty starburst galaxies, 
and AGN (at any redshift) would satisfy this extremely 
red criterion. As we shall demonstrate in the remainder 
of this letter, this indeed turns out to be the case. Note 
that late-type CWW Scd (cyan) and irregular (purple) 
galaxies are never sufficiently red to be classified as an 
ERO at any redshift^ and one would not expect to find 
any CWW late- type Sbc (green) EROs at 2; < 1.5. 

4. REDSHIFT DISTRIBUTION OF THE SPECTROSCOPIC 
ERO SAMPLE 

There are several thousand EROs in the EGS field. 
Here, we carry out a pilot study of the 87 EROs with 
confirmed DEEP2 spectroscopic redshifts. We use spec- 
troscopy in this letter only for redshift determination. 

The solid black histogram in the far left panel of Fig- 
ure [2 shows the redshift distribution of the 87 EROs in 
our sample. The EROs occupy a relatively narrow red- 
shift range at 0.76 < z < 1.42. The dotted black his- 
togram shows the redshift distribution of those galaxies 
in the DEEP2 EGS field with good quality spectroscopic 
redshift determinations (scaled down by a factor of 15). 

The fact that none of the EROs are located at 2; < 0.76 
is not a DEEP2 selection effect; in EGS, DEEP2 samples 
the fuh redshift range < z < 1.45 (|Fabeill2006HWillmei< 
I2OO6). It is caused by the additional ERO R- [3.6] > 4.0 
selection requirement, which effectively excludes low red- 
shift galaxies (Figured). The DEEP2+ER0 selection 
function also limits this sample to a fairly bright but nar- 
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row range of 18.5 < [3.6] < 20, with the one ex ceptio n 
of 12007954, a 17th magnitude AGN (|Le Floc^hll200l . 
Note that by selecting only those EROs with spectro- 
scopic redshifts, we may bias our sample against inclu- 
sion of any low-redshift extremely dusty galaxies (since 
optically faint galaxies will not pass the R < 24.1 selec- 
tion requirement of the DEEP2 survey). 

5. ERO CLASSIFICATION USING [3.6]-[8.0] COLOR 

It is possible to divide the EROs in this sample into 
three broad classes by means of their IRAC [3.6]- [8.0] 
color. In this section, we introduce each of these classes, 
show examples of their SEDs and discuss the sometime 
subtle features which are only apparent from such broad- 
baseline datasets. We also study their morphologies us- 
ing ACS images. 

From Figure^ early-type galaxies are expected to have 
the bluest ([3.6] — [8.0] —1.0) colors, dusty starburst 
galaxies intermediate colors ([3.6] — [8.0] o:^ —0.5), and 
AGN reddest colors ([3.6] - [8.0] +1.0). The right 
panel of Figure [3 corresponds to the dashed rectangle 
in Figure Q and shows the subsample of 87 EROs with 
spectroscopic redshifts. We classify [3.6] — [8.0] < —0.75 
galaxies as bulge-dominated early types (blue circles), 
-0.75 < [3.6] - [8.0] > 0.0 galaxies as dusty starburst^^ 
(red triangles) and [3.6] — [8.0] < 0.0 galaxies as power- 
law AGN types (orange boxes). As we shall demonstrate, 
these two simple color cuts do separate the classes rather 
well. This is because the redshift distribution of this par- 
ticular sample is rather narrow. For a sample spanning a 
broader redshift distribution, one might imagine utilizing 
a more complex color selection to isolate the meanderings 
of the different populations through color-color space in 
Figure ^ 

We classify 53 of the EROs as early types. They are 
located at 0.82 < z < 1.28 and their redshift distribu- 
tion is shown by the blue hatching in the middle panel 
of Figure [21 We classify 32 EROs as being dusty star- 
burst types. The dusty starburst population (shown by 
the red hatching) has a slightly wider redshift distribu- 
tion from 0.76 < z < 1.42 (presumably because galax- 
ies in this class can contain arbitrary amounts of dust). 
We classify two EROs as being predominantly power- 
law types. These are 12007878 at z = 0.99 and 12007954 
at z = 1.15, a nd are shown in solid orange (see also 
lKonidarisl[200l . 

6. EXAMPLE SPECTRAL ENERGY DISTRIBUTIONS AND 
ACS IMAGES 

Figures 01 and El (Plate 1) show examples of observed 
SEDs and 3'' x 3'' /-band (F814W) postage stamp ACS 
images. For comparison, also shown are CWW E (blue) 
and M82 (red) templates (§ O as they would appear at 
the redshift of each ERO. At z ~ 1, IRAC's 8/im chan- 
nel measures any emission from small hot dust grains 
(indicative of a dusty starburst galaxy) while the 3.6/im 
channel measures the stellar peak. A relative red (or 
blue) [3.6] — [8.0] color, therefore, can be used to dis- 
criminate for (or against) dusty starburst galaxies (and 
AGN) at this redshift. 

We use the word "starburst" to mean a galaxy whose IR 
luminosity is powered by star formation. We do not intend to imply 
these galaxies have total IR luminosity L(8 — lOOO/xm) < 10^^ Lq. 
Indeed, as we shall discuss in § Q most of these dusty starburst 
EROs would actually be classified as LIRGs. 
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6.1. Early-Type EROs 

Figureini (Plate 1) shows six examples of galaxies IRAC 
color-classified as early-type EROs. The CWW E tem- 
plate approximates the 13048898, 13019047, 13004276, 
and 12004426 SEDs very weh. The CWW E template 
appears to underestimate the SEDs of 13019309 and 
12008091 in the UV. In ah cases the ACS images show 
clear evidence of bulge-dominated morphologies, sup- 
porting our color classification as early types. 

A total of 35 of the 53 early-type EROs fall within the 
ACS footprint. A close inspection of all of the available 
ACS images reveals them to be predominantly bulge- 
dominated (E, SO or Sa-type spirals). A small number 
(~ 10%) appear to be undergoing mergers. We conclude 
that choosing a blue [3.6]- [8.0] color successfully selects 
for z = 1 EROs with old stellar populations and against 
those with dusty starburst or power-law features. 

6.2. Dusty Starburst and Power-Law EROs 

Figure ^ (Plate 1) shows five examples of galaxies 
IRAC color-classified as dusty starburst EROs, and one 
power-law ERO (12007878). The M82 template approx- 
imates the SEDs of sources 13042940, 12100899, and 
12007831 very weh, and that of 13004291 reasonably weh. 
In addition to being more mid-IR luminous than an early- 
type ERO of similar 3.6/im magnitude, the SEDs of these 
four sources are more luminous in the optical. Their ACS 
images show clear evidence of disturbed, peculiar, inter- 
acting or merging galaxies. 

Neither template approximates well the SED of 
12008048. Close inspection of the ACS image reveals 
it to be a face-on spiral galaxy. Since a normal spiral at 
z = 0.9 would not meet the R— [3.6] redness criterion, we 
infer that 12008048 must be an especially dusty spiral. 

A total of 17 of the 32 EROs we classify as dusty star- 
bursts fall within the ACS footprint. Close inspection 
of all of their images reveals many of them (~ 60%) to 
be disturbed or interacting galaxies, and the remainder 
(^ 40%) to be late-type spirals. 

12007878 shows a monotonically increasing power-law 
SED. Its host galaxy is clearly bulge-dominated and it 
can be unequivocally classified as an AGN-dominated 
source. 

7. DISCUSSION 

Turning to those EROs which have 24/im detections, 
there are five (17.4 < [24] < 19.1) amongst the 53 early- 
type and 29 (16.8 < [24] < 19.1) amongst the 32 dusty 
starburst EROs. Both power-law EROs have a 24/im de- 
tection (14.5 and 16.5). In four out of the five cases of 
early-type EROs with a 24/im detection, an ACS image 
is available. In all four cases, a merging galaxy with sev- 
eral close neighbors is apparent. The MIPS point-spread 
function is 6^' so it is possible that the signal in each case 
might be associated wit h another galaxy entirely. 

We use lLe Floc'h et all ((200 5. Fig. 7) to translate our 
24/im fiux limit (80uJy, 19.1 AB, 5cr) into a lower limit on 
the total detectable IR luminosity of an ERO in our red- 
shift range. Depending on the model used, at z ~ 0.75, 
we are capable of detecting sources more luminous than 
^ 5 X lO^^L© i.e., some dusty starbursts, and all LIRGS 
and ULIRGS. At 2; ~ 1.4, we are capable of detecting 
sources more luminous than ~ 5 x lO^^L© i.e., LIRGS 



and ULIRGS. Based on their actual 24/im magnitudes 
and redshifts, most of our "starburst" sources have in- 
trinsic luminosities between 10^^ and 10^^ L©, i.e., LIRG 
class. Both of the "AGN" sources have intrinsic luminosi- 
ties of > lO^^L©, putting them in the ULIRG class. 

Although, in this study, we utilized optical/IRAC col- 
ors to classify our ERO sample, IRAC- [24] color selection 
may be a m ore effective discriminator a t higher redshi ft 
(|Lacv et al . 2004; Ivison et al. 2004; Saii na et al.ll200.^ . 

Putting IRAC ERO colors in the context of perhaps 
more familiar rest frame optical colors. Figure shows 
restframe U — B color versus observed [3.6] — [8.0] color 
for all IRAC-detected galaxies with good [8.0] and good 
CFHT BRI photometry (see Willmer et al. 2006 for de- 
tails of the conversion to restframe U — B). The symbols 
are as in Figure [2l The filled circles in Figure [3 denote 
the 36 EROs which are also detected at 24/im (see also 
Figure [2j. The EROs we classify as early types gener- 
ally tend to lie on the U — B defined "red-sequence" and 
the EROs we classify as dusty starburst tend to lie in 
t^ "green valleyj'^ and at the top of the "blue cloud" 

8. CONCLUSIONS 

In this paper, we utilized the AEGIS dataset to ex- 
plore the nature of 87 EROs with confirmed spectro- 
scopic redshifts. IRAC imaging facilitated dividing this 
sample into three distinct classes using a simple [3.6] - 
[8.0] color selection technique. We presented SEDs and 
high-resolution ACS images supporting our (early-type, 
dusty starburst or power-law) classification. 

We showed that the three classes of ERO and their 
redshift distribution were broadly consistent with what 
one would expect from a simple color-color diagram 
(Figure We found 53 early- types, 32 dusty star- 
bursts/LIRGs, and two obvious AGN. Both of these 
AGN would be classified as ULIRGs. 

The selection of this particular sample of EROs was 
subject to the spectroscopic biases discussed in § ^ 
We next plan to extend this study to the several thou- 
sand IRAC-selected EROs in the AEGIS without spec- 
troscopic redshifts. 

This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. Support for this 
work was provided by NASA through an award issued 
by JPL/Caltech. ALC is supported by NASA through 
Hubble Fellowship grant HF-01182.01-A. We thank the 
referee. Matt Malkan, for constructive comments. GW 
thanks Mark Lacy, Adam Muzzin, Jason Surace, Ian 
Small and Mike Hudson for useful discussions, and UC 
Berkeley and UC Santa Cruz for their hospitality. 
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Fig. 1. — R — [3.6] versus [3.6] — [8.0] color-color diagram for IRAC-selected galaxies. The solid lines show non-evolving templates in 
color-color space for CWW early- type (blue), CWW Sbc (green), CWW Scd (cyan), CWW Im (purple), dusty starburst M82 (red), dusty 
starburst/ULIRG Arp220 (yellow), NGC 1068 (pink) and NGC 5506 (orange). The redshift is indicated for each type. A distinctive swath 
of galaxies is apparent curving from the left to the top of the Figure. The CWW early-type (blue) track follows this extremely well to 
z ~ 0.7, when the template noticeably begins to diverge from the "bluer" data. The plume in the lower left corner is caused by stellar 
contamination. Galaxies with R — [3.6] > 4.0 (i.e., above the black line) are defined as EROs. One would expect only high redshift (z > 0.8) 
early-type, (z > 1.0) dusty starburst galaxies, (z > 1.5) CWW Sbc, and AGN (at any redshift) to satisfy this extremely red criterion. The 
dashed black rectangle indicates the region of color-color space populated by the 87 EROs with confirmed spectroscopic redshifts (see also 
right panel of Figure |21 and §[3. 
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Z Z [3.6] - [8.0] 

Fig. 2. — The solid black histogram in the far left panel shows the ERO redshift distribution N(z) (which peaks for this sample of 87 
EROs at z = 1.15). The dotted black line illustrates the underlying DEEP2 selection function in the EGS field before the ERO criterion is 
applied. The dotted black line shows the redshift distribution N(z) of all 6023 galaxies with good quality redshift determinations, scaled 
down by a factor of 15 for comparison with the ERO N(z). The right panel (dashed black rectangle in Figure^} shows our IRAC [3.6] — [8.0] 
color division into early type (blue circles), dusty starburst (red triangles) and power-law (orange boxes) EROs. Filled symbols indicate 
the 36 EROs with a 24/im detection (19.1 AB, 5cr, §[7|. Also shown in the middle panel (and discussed in §[3, is the redshift distribution 
of the full sample (black) and the three sub-classes (using the same blue, red and orange color scheme). 
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Fig. 3. — Examples of observed SEDs (left) and ?>" x ?i" /-band ACS images (right) for six IRAC color-selected early- type EROs. A 
maximum of 12 passbands (^i^ B, g' , Subaru R, CFHT R, /, K [3.6], [4.5] [5.8], [8.0] and [24.0]) are shown. The tip of the downward 
pointing arrow indicates upper limits in the case of non-detections. Each panel shows the DEEP2 ID number, and the spectroscopic redshift. 
Also shown are a CWW E (blue) and M82 (red) template (§E| as they would appear at the redshift of each ERO. See § [H] discussion. 
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Fig. 4. As for Figure|3l but for five IRAC color-selected dusty starburst EROs and one power-law ERO (12007878). The M82 dusty 
starburst template approximates the 13042940, 12100899, and 12007831 SEDs very well, and that of 13004291 reasonably well. In most 
cases the ACS images show clear evidence of interactions and mergers. 12008048 is a face-on dusty spiral galaxy. 
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Observed [3.6] - [8.0] 

Fig. 5. — Restframe U — B color versus observed [3.6] — [8.0] color for all IRAC-detected galaxies with good [8.0] and good CFHT BRI 
photometry. Note: {U — B)vega = {U — B)ab - 0.85. Blue circles indicate bulge-dominated early-type EROs. Red triangles indicate dusty 
starburst EROs. Orange boxes indicate the two AGN. Filled symbols indicate a 24/im detection. Is is possible that the 24/im detections 
associated with each of the early- type EROs may be spurious, and actually associated with a close neighbor galaxy (see §|7|. The early- type 
EROs generally tend to lie on the U — B defined red-sequence and the dusty starburst EROs lie in the "green valley" and at the top of the 
"blue cloud". 



